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ABSTRACT 

Taking advantage of a high-resolution simulation coupled with a state-of-art semi- 
analytic model of galaxy formation, we probe the mass segregation of galaxies in 
groups and clusters, focusing on which physical mechanisms are driving it. We find 
evidence of mass segregation in groups and clusters up to the virial radius, both looking 
at the galaxy stellar mass and subhalo mass. The physical mechanism responsible 
for that is consistent with dynamical friction, a drag-force that brings more massive 
galaxies faster towards the innermost regions of the halo. At odds with observational 
results, we do not find the inclusion of low-mass galaxies in the samples, down to 
stellar mass M* = 10 9 Mq, to change the overall trend shown by intermediate and 
massive galaxies. Moreover, stellar stripping as well as the growth of galaxies after 
their accretion, do not contribute either in shaping mass segregation or mixing the 
radial mass distribution. Beyond the virial radius we find an “anti-mass segregation” in 
groups that progressively weakens in clusters. The continuous accretion of new objects 
and recent merger events play a different role depending on the halo mass onto which 
accreting material is falling. 
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1 INTRODUCTION 


The current scenario of galaxy formation can be divided in 
two distinct galaxy evolutionary families: the nature sce¬ 
nario, where galaxy properties depend on the state of the 
galaxies at the time of their formation, and the nurture sce¬ 
nario, where the galaxy properties are instead the product 
of environmentally driven processes that take place after the 
galaxy is accreted on la rge systems, such as g roups and clus¬ 
ters. As highlighted by !Pe Lucia et al.l (120121 1 . galaxies ex¬ 
perience several different kinds of environment during their 
lifetime, thus leading to the conclusion that the aforemen¬ 
tioned scenarios are twisted together. The mass distribution 
of galaxies in different environments is therefore the result 
of many physical processes that operate during galaxy evo¬ 
lution, such as gas cooling, star formation, stripping and 
mergers, and the environment is bel ieved to have a key-role 
in shaping galaxy properties (e.g.. 
Ta^ka^t^aJjgOO^ jBEntmyet al.1 


Kauffmann et all 200J; 

20051 : IWeinmann et al l 


20061 : Blanton fe Berlindf 2007 ). We now know that galaxy 


properties such as mass, colour, gas content and age are 
strongly related to the environment, such that galaxies in 
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denser environments are typically more massive, redder, less 
gas-rich and older. 

Dynamical friction is a key factor in the link between 
galaxy growth and environment. This has been pointed- 
out by several aut hors in the context of dark matter sub - 
structures (see e.g. iDe Lucia et al.l[2004l : IContini et al.l[2012h 
as well as the connec tion between galaxy stellar mass 
and subhalo ma ss (e.g ., Ipao et "all 12004 IWang et al.l 120061 : 
lYale fe Ostrikerl l2006h ■ According to the dynamical fric¬ 
tion formula (given bv lChandrasekhanll943ll . more massive 
galaxies fall faster towards the innermost region of the clus¬ 
ter, thus leading to a segregation in mass. Under this pic¬ 
ture, the characteristic gala xy mass is expected to increase 
with increasing halo mass IlYang et ahl 120051 : IZheng et alJ 
120051 : Ivan den Bosch et aLll2008i f because galaxies are more 
massive in larger objects, and more massive galaxies should 
preferentially reside in the centre. Hence, if dynamical fric¬ 
tion plays an important role, its effect on the distribution of 
galaxies in large system must be observable. On the other 
hand, stellar stripping due to tidal forces between galaxies 
and the cluster potential, as well as ongoing star formation 
in galaxies as they fall towards the center, might affect the 
distribution in mass. Tidal stripping removes part of the 
stellar mass of galax ies, and is expected to be more effective 
in massive galaxies jRudick et al.l 120091 : iMartel et al.l 120121 : 
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Watson fe Conrovl 120131 : lLaporte et al.l 120131 : IContini et al.l 
2014), while star formation leads to an increase in stellar 
mass. 

Despite the large amount of observations and effort 
spent on the topic, there is not yet a common agreement 
in the literature. There i s evidence in favour of mass segre¬ 
gatio n in large systems lLa res_eLal] HixMl; |McIntoshet_alJ 


20051: van den Bosch et al.l 12 0081; Presotto et al.1 20121 : 
Baloeh et, al.1120141 : iR.oberts et al.l 1201 fil l. but at the same 
time there is signi ficant evidence for a lack of mass segre¬ 
gation in clusters (von_der Linden et al.ll2Qlol : IVulcani et al.l 


120131 : IZioaro et al. 1201 .fl ). Thus, the controversy over mass 
segregation re mains unsol ved._ 

Recently. iRoberts et al.1 d2015l l. using group catalogues 
derived from th e Sloan Digit al Sky Survey Data Release 7 
(SDSS DR7. lAbazaiian et al.l [20091 1. find evidence for mass 
segregation out to 2 -i? 200 , with strength that scales inversely 
with halo mass. These authors argue that the conflicting 
results in the literature might be reconciled by considering 
that mass segregation appears more evident once low-mass 
galaxies are included in the sample, and that it is a function 
of halo mass, with clusters showing little or no segregation. 
This points-out the need to probe a large range of halo mass 
and different thresholds in galaxy stellar mass, and the need 
of timely investigations using theoretical models. 

Theoretical models have n ot suffic iently c ontrib uted to 
the debate. Only very recently IVulcani et al.l d2014l l (here¬ 
after V1 4), by means of semi-anal ytic models of g alaxy for¬ 
mation m. Lucia fe Blaizoti 120071 and IGuo et ai]|201ll 'l. in¬ 
vestigate the galaxy stellar mass function as a function of 
the environment. Although they do not explicitly focus on 
mass segregation, these authors find that the stellar mass 
function does not strongly depend on the cluster-centric dis¬ 
tance, although there is a hint of mass segregation in low and 
intermediate halo masses. 

In this letter we focus on the role played by dynami¬ 
cal friction on the observed mass segregation in groups and 
clusters. We take advantage of a state-of-art semi-analytic 
model which includes a refined treatment of stellar strip¬ 
ping that can be switched-on/off. Our simulation allows to 
achieve the goals highlighted above, since we can probe a 
large range of halo mass and use different thresholds in stel¬ 
lar mass in order to: a) investigate the dependence of mass 
segregation with halo mass; b) understand if it is driven by 
the inclusion of low-mass galaxies in the sample; c) isolate 
the role of stellar stripping; d) isolate the growth of galaxies 
after accretion. 


2 METHODS 

The simulation used in this pa per an d the semi-analytical 

model are based on iKang et al.1 (|2012l l . We refer the read¬ 

ers to that paper for details. Here we simply introduce the 

main prescripti ons. The simu lation was performed using 

Gadget-2 code ISpr ingel 2005) with cosmo logical parame¬ 

ters a dopted from the WMAP7 data release (IKomatsu et all 
l201lT) . namely: Q x = 0.73, = 0.27, Q b = 0.044, h = 

0.7 and os = 0.81. The simulation box is 200 Mpc/h 
on each side using 1024 2 3 * * * * particles, each with mass 5.64 • 

10 8 Mq/U 1 . The merger trees are constructed by follow¬ 

ing the subhaloes resolved in FOF haloes at each snap- 


Table 1. Our simulated haloes have been split in four samples, 
according to the halo mass. In the first column, we give the name 
of the sample, the second column indicates the range of M 200 
values corresponding to each sample, while the third gives the 
number of haloes in each sample. 


Sample 

Mass Range (M 200 ) 

Numb, of Haloes 

Small Groups 

[10 13 - 5 • 10 13 ]A7 o /i _1 

2286 

Large Groups 

[5 ■ 10 13 - 10 14 ] A7 0 h —1 

275 

Small Clusters 

[10 14 - 5 ■ 10 14 ]M 0 /i _1 

146 

Large Clusters 

> 5 ■ 10 14 Mq1i _1 

9 


shot (e.g., IKang et al.l 2005h m aking use of the algorithm 
SUBFIND ( Soringel et al.1200ll f. The semi-analytic model is 
then grafted on the merger trees and self-consistently mod¬ 
els the physical processes governing galaxy formation, such 
as gas cooling, star formation, supernova and active galac¬ 
tic nucleus feedback. The galaxy luminosity and colours are 
ca lculated based on the ste llar population synthesis model 
of lBruzual fc Gharlotl d2003h u sing a Chabrier stellar initial 
mass function I Chabrierl l2003l. Although we select galaxies 
with stellar mass larger than 10 9 Mq/i _1 , it is worth not¬ 
ing that the stellar mass completeness of the galaxy cata¬ 
logue provided by the model run on the simulation is about 
10 8 * Moh- 1 . 


We updated the semi-analytic model described above 
by adding t he prescription Mod el Tidal Radius+Mergers 
presented in IContini et al.l (12014I - ). which accounts for the 
formation of the intra-cluster light via stripping processes 
and mergers. The stellar stripping channel accounts for tidal 
forces that might strip a given fraction of stellar mass from 
satellites galaxies as they approach the innermost regions 
of dense environments such as galaxy clusters. Depending 
on the strength of the tidal force, the galaxy might be to¬ 
tally disrupted. The merger channel, instead, considers that 
a fraction (20 per cent in that model) of the satellite stellar 
mass gets unbound in violent relaxation processes t hat hap - 
pen during galaxy-galaxy mergers. In IContini et al.l (120141 1 
we have verified that such a simple prescriptions reproduces 
approximately the res ults of the numerical simulations by 
IVillalobos et al.1 |2012l l. It is however worth noting that, in 
reality, the fraction of star s that are unbound sho uld depend 
on the orbital circularity dVillalobos et al.ll2014l l. 


3 RESULTS 

In order to examine environmental dependences we study 

4 samples of haloes, ranging from small groups (with mass 
M 200 in the range [10 13 — 5 • 10 13 ]Mq/i - 1 ), to large clus¬ 
ters (with mass larger than 5 • 10 14 A7q/i _1 ). Details of the 
samples are given in Table [T] 

In the top panels of Figure [T| we plot the mean stellar 
mass of galaxies as a function of radial distance considering 
the mass at redshift z = 0 (solid lines) and at the time of 
accretion Q (dotted lines) for our 4 samples of haloes (dif¬ 
ferent colours), considering all satellite galaxies (explicitly 


1 The time of accretion of a given galaxy is defined as the last 
time the galaxy is central. 
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Figure 1. Mean stellar mass of galaxies as a function of radial distance considering the mass at redshift 2 = 0 (solid lines) and at the 
time of accretion (dotted lines). Dashed lines represent la statistical errors computed for the mass at redshift 2 = 0 (we omit the scatter 
around the mean mass at accretion for a more readable plot). All galaxies having stellar mass larger than 10 10 Mq belonging to the 
FOF-group are considered (top panels), for halo masses within [10 13 — 5 • 10 13 ]Mgh _1 (blue lines), [5 • 10 13 — 10 14 1 M^fe^ 1 (red lines) , 
[10 14 — 5 • IQ 14 ] JVf © h ~ 1 ( green lines), and masses > 5 • 10 14 M©/i —1 (purple lines). Symbols represent data points of lVulcani et al]|2014l . 
which are predictions of lDe Lucia fc BIaizonl2007l model. The bottom panels show the same information for the same samples of haloes, 
but considering all galaxies within a sphere of radius 2 ■ R 200 centred on the halo centre. 




Figure 2. Mean dark matter mass in subhaloes of galaxies with mass larger than IO 10 M©h ~ 1 in the same samples of haloes shown in 
the bottom panels of Figure |T] Due to the presence of ’’orphan” galaxies in our model (see text), this population of galaxies are missed 
here. 


excluding centrals) with stellar mass larger than 10 10 Mg 
belonging to the FOF-group. We find a clear mass segre¬ 
gation up to the virial radius, with more massive galaxies 
preferentially located in the innermost regions, close to the 
centre. To better quantify the strength of the segregation we 
compute a linear fit out to the virial radius for all samples 
and show the result in Table [2] The zero-point ( b ) decreases 
from small groups to big clusters of about 0.2 dex, while the 


slope (a) flattens significantly from -0.678 in small groups 
to -0.227 in big clusters. Moreover, none of the slopes are 
consistent with zero even within 3 ■ a. 

Interestingly, at r ~ R200 we find an upturn, and the 
mean stellar mass starts to increase again up to 2 • R200, 
from small groups to small clusters. As shown by the plot, 
the result still holds when considering the stellar mass at 
the time of accretion. In this case the trend is the same, 
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Table 2. Slopes, zero-points and their errors obtained with a 
linear fit, log M* = a ■ ( R/R 200 ) 4- 6, done out to the virial radius, 
for all samples shown in the upper panels of Figure [Tj 


Sample 

a 

<y a 

b 

<*b 

Small Groups 

-0.678 

0.012 

10.797 

0.005 

Large Groups 

-0.483 

0.023 

10.728 

0.011 

Small Clusters 

-0.339 

0.021 

10.665 

0.011 

Large Clusters 

-0.227 

0.045 

10.681 

0.023 


but the mean stellar mass at each radial distance is slightly 
lower. This means that galaxies have had time to grow in 
mass by redshift z = 0, but also that the effects of stellar 
stripping and star formation do not drive trends in mass 
segregation, suggesting that the main driver is of dynamical 
na ture. We compare our res ults with results of V14 given 
bv lDe Lucia fe Blaizotll2007l ’s model (black symbols in the 
top panels of Figure [lj. Despite these points are shifted-low 
with respect to ours, they clearly show the same trend and 
a hint of upturn in the external regions in small and large 
groups, while they suggest neither mass segregation within 
the virial radius nor upturn in large clusters. 

In the bottom panels of Figure Q] we show the same 
information plotted in the top panels, now considering all 
galaxies within a sphere centred in the centre of the halo 
and having a radius equal to 2 • R200 ■ This procedure results 
in a sample including all of those galaxies in the top panels of 
Figure [T] as well as all galaxies within the sphere that are not 
gravitationally bound to the group/cluster, and results in a 
sample that is more comparable to observational samples. 
A comparison between the top and bottom panels indicates 
that adding galaxies not belonging to the FOF-group, but 
within 2 • R200, attenuates the steepness of the stellar mass- 
radial distance relation beyond R200 in groups, and basically 
flattens it in clusters. 

In Figure[2]we show the dark matter mass (DM) in sub¬ 
haloes associated with galaxies with stellar mass larger than 
10 10 M 0 r', for our samples of haloes. In our semi-analytic 
model we keep following the so-called ’’orphan-galaxies”, i.e. 
those galaxies that have lost their subhalo after accretion in 
a larger system. In Figure [2] we take the samples of galaxies 
considered in the bottom panels of Figure |T] without the 
orphan-galaxies. The trend of mass with radius shown by 
DM is the same found for stellar mass, a clear mass segrega¬ 
tion up to the virial radius and an upturn beyond it, in all 
samples except for the small groups. This suggests that mass 
segregation might be driven mostly by dark matter rather 
than by physical processes involving only stellar mass. 

The prediction of our model is then consistent with ob¬ 
servations that find evidence for mas s segre g ation in groups 
and clusters ( van den Bosch et al.l 20081: Presotto et al.l 
120121 : iBaloeh et al.l 20141 : Roberts et al.l 2015T ). at least up 
to the virial radius. The upturn beyond the virial radius 
disappears on cluster scales, particularly when the sample 
includes non-FOF galaxies. We repeated the same analysis 
including low-mass galaxies (down to 10 9 Mq). The inclusion 
of low-mass galaxies only changes the mean stellar mass at 
each radial distance, but does not significantly change the 
overall trend. 


4 DISCUSSION 

Mass segregation in evolved objects is generally believed to 
be a consequence of dynamical friction, which causes more 
massive galaxies to fall to the centre of haloes more quickly. 
Nevertheless, the assembly history of haloes, as well as other 
physical processes taking place in dense environments, might 
play an important role. Galaxies can grow in stellar mass via 
star formation or mergers, and the rate of the growth could 
depend on their stellar mass. In addition, tidal forces in 
groups and clusters cause stellar stripping. Although dynam¬ 
ical friction brings massive galaxies towards the centre, this 
is also t h e envi ronment in which tidal stripping gets stronger 
(jContini et al.ll2014ll . It is possible that tidal stripping could 
be strong enough to mix the population of galaxies and no 
mass segregation would be present. We investigated on this 
switching-off our stellar stripping and merger channels for 
the formation of the intra-cluster light, and found similar re¬ 
sults as those reported in Figure Q] (plots not shown), mean¬ 
ing that stellar stripping is not enough for mixing the galaxy 
distributio n. This is confirmed by V 14’s data, that are pre¬ 
dictions of lDe Lucia fe Blaizotll2007l models, which does not 
have any prescription for stellar stripping. 

As seen in Section [3] the top panels of Figure [T] show 
that beyond -R 200 the stellar mass-radial distance relation 
changes slope until it flattens in clusters. When considering 
all galaxies (bottom panels of Figure [lj in the sphere cen¬ 
tred on the centre of the halo, the slope flattens at lower halo 
mass, in large groups. This may be due to pre-processing and 
accretion of new galaxies by the FOF-group, which could 
reduce the strength of ’’anti-mass segregation” found out¬ 
side the virial radius, that progressively weakens from small 
groups to large clusters. This is supported by the result 
found in Figure [5] for DM subhaloes. Accretion of smaller 
haloes and mergers of objects with comparable mass dur¬ 
ing the assembly of each halo have different consequences 
depending on the virial mass of the object that they are in¬ 
falling onto or merging with. Large groups and clusters form 
later, and this may lead to the upturn due to recent mergers 
of DM haloes. In terms of stellar mass we do not find an up¬ 
turn on cluster scales due to the non-linear relation between 
DM and stellar mass. This means that, in our model, the 
stellar mass growth in massive haloes is not as steep as the 
DM mass growth, and that the star formation efficiency is 
higher in low-mass haloes with respect to larger ones. This 
is consistent with the predicted stellar mass function, from 
which it is clear that semi-analytic m odels find more low- 
mass galaxies than observed (see, e.g.. iGuo et al. ism and 
references therein). 

We investigate the relationship between dynamical fric¬ 
tion and halo mass in Figure [3] where we plot the residual 
merging time at z = 0 versus the radial distance, for the 
same samples of haloes and considering all galaxies in the 
sphere centred on the centre of the halo and having a radius 
equal to 2 • R200, as done in the bottom panels of Figure [l] 
The residual merging time is simply derived as the dynam¬ 
ical friction timescale minus the time elapsed since accre¬ 
tion. The dynamical friction timescal e has b een evaluated 
at the time of accretion and following 1 Jiang et al.ll200j for¬ 
mula (eqn. 5 of their paper). Two interesting features arise 
from this plot. First, galaxies that are closer to the centre 
are also those with the smallest residual merging time and, 
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Figure 3. Residual merging time at z = 0 (derived as the dynamical friction timescale minus the time elapsed since accretion) as a 
function of their radial distance from the halo centre. All galaxies within a sphere centred on the halo centre and having a radius equal 
to 2 ■ f?200 are considered. 


overall, the relations fit well with the behaviour shown in 
the bottom panels of Figure |T] Second, the residual merging 
times increase from groups to clusters. This means that mass 
segregation in groups can be explained by invoking dynami¬ 
cal friction, which has enough time to bring massive galaxies 
in the innermost regions of the group. This still applies in 
clusters, where the typical dynamical friction timescale is 
long enough to distribute galaxies according to their mass, 
but on average much longer than in groups (even longer than 
a Hubble time). 

We are aware that the ’’anti-mass segreg ation” we find 
beyond the virial radius is not observed (e.g. iRoberts et ali 
120151) . There might be several reasons for that, including 
contamination of foreground and background galaxies. We 
do not have any contamination of such objects by construc¬ 
tion. On the other hand, when we include unbound galaxies 
as done in the bottom panels of Figure |T] instead of FOF 
galaxies only (as done in the top panels of the same fig¬ 
ure), the relation considerably flattens. As argued above, 
this might be due to recent mergers between DM haloes, 
which could reduce the strength of anti-mass segregation 
found outside the virial radius. The upturn beyond R 200 is 
remarkably sharp in groups, and appears sensitive to the 
inclusion of unbound galaxies. A full explanation of the up¬ 
turn deserves a more detailed analysis, which must take into 
account the merging histories of the haloes considered. We 
plan to address this point in a future paper. 


5 CONCLUSION 

In this letter we focus on the mass segregation in simulated 
groups and clusters coupled with a state-of-art semi-analytic 
model of galaxy formation that accounts for the physics of 
baryons. 

We find a non-negligible mass segregation in groups and 
clusters up to the virial radius, and the level of segregation 
is statistically significant at all halo masses. The strength of 
the segregation is found to be a function of halo mass, such 
that is weaker at higher halo mass. This can be explained 
by invoking dynamical friction, which brings more massive 
galaxies faster towards the innermost regions of the halo. 
The mass segregation trends are insensitive to both the in¬ 


clusion of low-mass galaxies in the sample (down to stellar 
mass M* = 10 9 Mg) and stellar stripping. 

Moreover, beyond the virial radius we find an ’’anti¬ 
mass segregation” in groups that progressively weakens in 
clusters. The upturn beyond the virial radius highlights that 
these galaxies are intrinsically more massive than those in 
the inner regions. Galaxies in the outskirts have been re¬ 
cently accreted and thus have been centrals for a longer time. 
They then had more time and more chance to grow in mass, 
due to cooling and enhanced star formation. Considering all 
galaxies in a sphere centred on the centre of the halo and 
having a radius equal to 2 ■ -R 200 (i.e. not only those belong¬ 
ing to the FOF-group), the trend beyond the virial radius 
weakens, approaching observational findings. 

Interestingly, mass segregation is found also by look¬ 
ing at the dark matter in subhaloes, with similar features 
shown by the stellar mass. This result, along with the fact 
that neither stellar stripping nor star formation (nor both 
together) shape mass segregation after accretion, points out 
that the main driver of the radial distribution of galaxies 
has a dynamical nature. We showed that dynamical friction 
is the most likely candidate, but, in order to quantify its rel¬ 
ative contribution in shaping mass segregation and confirm 
the nature of the upturns beyond the virial radius, a more 
detailed analysis is needed. 
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